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ABSTRACT 

We evaluate the achievable maximum energy of nuclei diffusively accelerated by shock wave in the 
jet of Cen A, based on an updated model involving the stochastic magnetic fields that are responsible 
for recent synchrotron X-ray measurements. For the maximum energy analysis, conceivable energy 
constraints from spatiotemporal scales are systematically considered for the jet- wide including discrete 
X-ray knots. We find that in the inner region within ~ 1 arcmin from galactic core, which includes 
knots AX and BX, proton and iron nucleus can be accelerated to 10 19 — 10 20 and 10 21 eV (10 — 100 EeV 
and ZeV) ranges, respectively. The upper cutoff energy of the very energetic neutrinos produced via 
photopion interaction is also provided. These are essential for identifying the acceleration site of the 
ultra-high energy cosmic ray detected in the Pierre Auger Observatory, which signifies the arrival from 
nearby galaxies including Cen A. 

Subject headings: acceleration of particles — galaxies: individual (Cen A) — galaxies: jets — magnetic 
fields — methods: analytical — turbulence 



1. INTRODUCTION 

To date, the 27 ultra-high energy cosmic ray (UHECR) 
events of the energy exceeding 5.7 x 10 19 eV have 
been detected in state-of-the-art Auger observatory; in 
particular, the anisotropy and significant correlation 
with the active galactic nuclei (AGNs) residing within 
75 Mpc, namely GZK horizon (Greisen 1966; Zatsepin 
& Kuzmin 1966), were discovered (Abraham et al. 
2007). Surprisingly, the results indicate that two of these 
UHECR events have arrived within 3° of CentaurusA 
(NGC5128), a closest galaxy (3.7 Mpc from us; Fer- 
rarese et al. 2007). Although the detailed position of the 
UHECR production site is still unresolved, the galactic 
core accompanied by a supermassive black hole, bipolar 
jets, giant radio lobes (Hardcastle et al. 2009), and so 
on, will be enumerated as the favored candidates. It is 
desired that genuine theoretical survey be expanded for 
clarifying the particle acceleration mechanism feasible at 
these sites, in order to provide the physical interpreta- 
tion of the observed intriguing results, which include the 
recent detection of high-energy gamma rays (Aharonian 
et al. 2009). 

When closely looking at the anatomy of the large-scale 
jets in the Cen A galaxy (see, Israel 1998, for a review), a 
clue to solve this challenging problem seems to be in the 
knotty regions resolved in a deep X-ray image (e.g., Kraft 
et al. 2002). In an inner region near the galactic core, 
the rugged features are associated with the shocks, which 
are considered to be formed where a flowing plasma col- 
lides with obstacles (Hardcastle et al. 2007). Accord- 
ing to the latest theory of plasma kinetic transport, it 
is known that such a colliding plasma is likely unstable 
for the electromagnetic current filamentation instability, 
which generates small-scale magnetic fluctuations with 
the order of plasma skin depth (Medvedev & Loeb 1999; 
Honda 2004 and references therein). In the nonlinear 
phase, the filaments preferentially coalesce one another, 
self-organizing larger scale filaments with stronger mag- 
netic fields. This reflects a stochastic nature of the in- 



verse energy cascade of plasma turbulence. The magnetic 
fluctuations are non-perturbative, and in the strong tur- 
bulence regime, in that the energy density is comparable 
to the thermal energy density of plasma bulk (Honda et 
al. 2000a,b). 

Indeed, recent high-resolution observations reveal that 
the jet of Cen A is in part dominated by the filamentary, 
sometimes edge-brightened features (Hardcastle et al. 
2007), as inspected in the well-confirmed jet of a nearby 
galaxy (Vir A; Owen et al. 1989). Worth noting thing is 
that inherent inner structure similar to these has been 
discovered in many other jets (e.g., Cyg A: Perley et al. 
1984; 3C 353: Swain et al. 1998; 3C 273: Lobanov & Zen- 
sus 2001; 3C438: Treichel et al. 2001; Mrk501: Piner et 
al. 2009). Also, it was recently found that a filamentary 
jet model naturally provides the comprehensive explana- 
tion for the complicated spectral variability observed in a 
TeV blazar object (Mrk421; Honda 2008). These ingre- 
dients strongly encourage us to take the inhomogeneous 
magnetic field effects into account for precisely modeling 
the Cen A jet as a cosmic-ray accelerator. 

For a standard, diffusive shock acceleration (DSA) sce- 
nario (e.g., Blandford & Eichler 1987), the strong tur- 
bulence plays an essential role in particle scatterers is 
expected. In the stochastic medium, the higher energy 
heavy particles tend to freely meander (Honda & Honda 
2005), albeit electrons (and positrons) are, if anything, 
likely gyro-trapped in the magnetized filaments, suffer- 
ing stronger radiative cooling (Section 3.1). From the 
fact that in the Cen A jet, electron synchrotron emis- 
sions appear to be, in part, already diffusive, it is thus 
inferred that for nuclei (say, proton) the conventional ap- 
proximation of small-angle resonant scattering will be in- 
adequate for describing the spatial diffusion, and instead, 
the three-dimensional rms deflection becomes rather fea- 
sible. Importantly, the latter facilitates the back and 
forth of particles across the shock, increasing the effi- 
ciency of DSA. Based on this notion, Honda & Honda 
(2004) have argued that a nucleus could be accelerated 
to the UHE range at a bright jet knot of Vir A, though the 
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arrival from the galaxy (~ 4 times more distant than the 
distance to CenA) is not yet signified (see, e.g., Stanev 
2008, for the useful discussions). 

Besides, the bound electrons are co-accelerated, to 
emit the synchrotron photons, which could serve as a tar- 
get of the accelerated protons. Then, the resonant inter- 
action of the UHE protons with the target photons could 
be a major neutrino production process via the decay 
channels triggered by photopionization: p-y — > tt^X — ► 
[i^v^ — ► e ± v e v ll (Romero et al. 1996). The emitted neu- 
trinos propagate on the straight, unlike the charged par- 
ticles that suffer, more or less, deflection by intergalac- 
tic magnetic fields (e.g., Vallee 2004), so that they play 
the role of a powerful and complementary messenger of 
the UHECR production, in light of the observability at 
a forthcoming kilometer neutrino telescope (e.g., Halzen 
& Hooper 2002). Recently, Cuoco & Hannestad (2008) 
have proposed the model spectra of high-energy neutri- 
nos from Cen A, but putting the detailed mechanisms of 
the particle accelerator into effect has remained unsolved. 

In this paper, based on the filament model, we esti- 
mate the maximum possible energies of a proton and 
iron diffusively accelerated in the Cen A jet, and also, the 
energy of the produced neutrinos, providing the encrgy- 
equipartition among the flavors. The conceivable mech- 
anisms of the energy restriction are taken into consid- 
eration for the jet wide including bright knots. It is 
addressed that the energy is limited dominantly by the 
shock operation time or particle escape loss, rather than 
radiative loss and nucleus-nucleus collision timescales. 
As a result, we find that in the inner region, which con- 
tains the X-ray knots BXn, proton and iron can be ac- 
celerated beyond the Auger limit. In particular, the ex- 
pected highest energy reaches, for iron, 3 ZeV and more 
(around knot BX5); that is to say, the CenA jet is a 
"Zevatron" worthy of the candidate (in addition to the 
Vir A, put forth by Honda & Honda (2004)). 

The present analysis virtually provides a substantially 
extended version of the previous simple analysis by 
Romero et al. (1996), and hence, a particular attention 
has been paid to highlight the new points including (1) 
the proposal of the improved theoretical model compati- 
ble with updated X-ray measurements (Section 2.1) and 
extended arguments on the turbulent magnetic field (Sec- 
tion 2.2), (2) the non-resonant diffusion scenario involved 
in the DSA (Section 3.1), and (3) thorough survey of the 
temporal (Section 3.2) and spatial (Section 3.3) limits, in 
order to figure out the maximum energies of proton and 
iron achievable in the inner region (Section 3.4), and es- 
timate neutrino energy (Section 3.5). The application of 
the maximum energy analysis to the outer region of the 
large-scale jet is also provided (Section 4). At last, the 
discussion on the comparison with the previous relevant 
results is expanded (Section 5). 

2. JET WIDE AS A CANDIDATE UHECR FACTORY 

2.1. The Kinematic and Spectral Properties 

Cen A is a giant elliptical galaxy (classified as Fanaroff- 
Riley (FR) type-I), which contains the bipolar jets 
ejected from the galactic core. The brighter jet extends 
toward the northeast direction beyond the projection of 
6* pro j ~ 250" (# P roj = 1" corresponds to 17 pc; Hardcastle 
et al. 2007), with its inclination with respect to the line 



of sight (the angle is referred to as the viewing angle, 
throughout). The apparent speed of motion is sublu- 
minal, exhibiting around 0.5c, where c is the speed of 
light (Hardcastle et al. 2003, 2006; Horiuchi et al. 2006; 
Brookes et al. 2006). Concerning the knots distributing 
along the jet, a model in which they are simply com- 
pressions in the fluid flow has been ruled out; and in- 
stead, they are considered to be privileged sites for the 
in situ particle acceleration (see, Hardcastle et al. 2003, 
for a detailed discussion). Although the velocity distri- 
bution inside the knot-like features is still not revealed, 
the knots may trace the stationary shocks, which could 
be formed as a result of interaction between the jet fluid 
and internal obstacles (Blandford & Konigl 1979a, b). For 
updated spectral indices around a r ~ 0.7 in observed ra- 
dio and radio-infrared continuum (Brookes et al. 2006), 
we read off the energy spectral index of energetic elec- 
trons of p = 2a r + 1 = 2.4. The radiation spectrum 
extending to the X-ray band is considered to be of the 
electron-synchrotron originator, which contains a key in- 
formation to understand the real mechanism of the accel- 
erator that simultaneously operates for ions, as long as 
ion abundance is finite in the jet interior (e.g., Evans et 
al. 2004; Markowitz et al. 2007). In this sense, below we 
deal with somewhat more details of the X-ray spectrum. 

Recent Chandra observations have successfully re- 
solved the X-ray knots smaller than the width of 
the radio-emitting column with diffuse X-ray emissions 
(Kraft et al. 2002; Hardcastle et al. 2007), although the 
inner structure of the knots remains unresolved as yet. 
When the filamentary morphology inspected in the jet 
(introduced in Section 1) retains a self-similar character- 
istic (Honda 2008), the magnification of the spatial scale 
of the X-ray knots could repetitiously reveal the compo- 
sition of smaller scale compactors and diffuse emitters. 
If this insight is correct, the DSA property is incurred by 
the filamentary turbulent state (Honda & Honda 2004, 
2005). Considering a similar situation that allows for 
magnetized filamentation, Spitkovsky (2008) have veri- 
fied that the energy gains truly occur as particles bounce 
between the upstream and downstream regions of a colli- 
sionless shock, and found p = 2.4 ± 0.1, amenable to the 
aforementioned value. 

According to the argument expanded by Honda & 
Honda (2007), the dominant synchrotron emissions 
(around the spectral break of <~ 10 14 Hz) are from the 
major electrons bound to the magnetized filaments with 
the maximum transverse size (i.e., the maximum coher- 
ence length). The electrons in smaller scale filaments 
down to a certain characteristic scale A c can be acceler- 
ated to higher energy, because of the weaker synchrotron 
loss in the weaker magnetic fields. In fine filaments 
smaller than A c , electron acceleration is limited by the es- 
cape from the filaments. 1 It is pointed out that the free 
electrons meandering in the magnetized filaments con- 
tribute to engender diffuse synchrotron photons (in the 
fundamental process). Anyhow, it follows that the spa- 
tial size scales of filaments are identified with the local 
maximum electron energies. In this context, particularly, 

Note here that this elementary process does not mean the es- 
cape (in meaning often said) from an entire system for particle 
confinement, like a jet column. The latter is considered in Sec- 
tion 3.3. 
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TABLE 1 

Summary of the Parameter Values 



Region 


^proi (P c ) a 


i?(pc) a 


Rout(pc) a 




/3 C 


AX1 


260 


352 


819 


260 


5/3 


AX2 


310 


352 


819 


250 


5/3 


AX3 


430 


352 


819 


290 


5/3 


AX4 


480 


352 


819 


230 


5/3 


BX2 


980 


819 


1310 


160 


5/3 


BX5 


1100 


819 


1310 


290 


5/3 



a cf. Kraft et al. (2002). 

b Equipartition values excluding the beaming effects (cf. 
tion 2.2). 

c The spectral index of the filamentary turbulence (cf. 
tion 2.1). 
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it appears that the electron synchrotron emissions from 
the major part, larger than A c , dominantly contribute 
to constitute an extended (typically X-ray) continuum 
of spectral energy distribution, instead of sharp cutoff. 
Taking into account of the flux density reduction due 
to the weak magnetic intensity, we can evaluate the F v 
spectral index above the break frequency, to give 



«x = 



(5/? -l)p-Q3 + 3) 
2(3/3+1) 



(1) 



where /3 corresponds to the filamentary turbulent spec- 
tral index. Here, we have assumed the Kolmogorov- 
type fluctuations superimposing on the local magnetic 
fields involved in the filaments. In Equation (1), we 
find that, for p = 2.4, the values of (3 = 5/3, 2, and 
3 (suggested by Honda (2008)) lead to a x = 1.1, 1.2, 
and 1.4, respectively. Interestingly, these are accommo- 
dated with the measured values of both a x = I.OO+q'^ 
in between knots A2 and B included in the projection of 
Voj < 60 " (Kataoka et al. 2006) and a x = 1.2 ± 0.2 in 
60" < 6» proj < 190" (Hardcastle et al. 2007). 

2.2. Turbulent Magnetic Fields 

It is expected that the magnetic field vectors are al- 
most randomly oriented inside the jet, and as is, the syn- 
chrotron radiation tends to be diffuse with the reduced 
polarization, with exception for that from the electrons 
bound to a pronounced coherent field. The field strength 
might be appropriately characterized by the rms value 
(denoted simply as B; Honda & Honda 2005), although 
it is hard to directly determine the exact value from inte- 
grating the turbulent spectral intensity over the jet wide. 
For simplicity, we here compare it with an equipartition 
value. Considering the emitting volume with the size of 
- 10" (Kraft et al. 2002; Hardcastle et al. 2006) that 
is larger than the size supposed by Kataoka & Stawarz 
(2005), we infer that the typical field strength may be 
represented by 100 )iG in the inner region including 
bright knots A and B, and less in the outer region (Hard- 
castle et al. 2006). In Table 1, the revised equipartition 
value, denoted as -Beq,5=i, for the simple case taking no 
rclativistic beaming effect into account is listed. The en- 
ergy density levels (u m ) are thought of as being larger 
than (or, at least comparable to) the radiation energy 
density (u rad ) (cf. Bai & Lee 2001). 

In the radial edge of the jet, uncanceled magnetic 
fields in the envelope of the clustered filaments coop- 
erate to apparently constitute a long magneto-tail ex- 




1 1 

Deprojected distance L (kpc) 



Fig. 1. — Opening angle of the Cen A jet (in the regions including 
the X-ray knots; labeled) <fi vs. the deprojected distance from 
galactic core L. Here, the definition of <fi = 2tan _1 (i? sin 0/L pro j ) 
has been used, where R and L pI0 j are the radius of jet and 
projected distance from the core, respectively (cf. Table 3), and 8 
is the viewing angle, assumed to be 15° (Hardcastle et al. 2003). 
The solid line indicates the interpolation among the points. 



tending to the radial direction (Honda & Honda 2004). 
The well-pronounced tail would have an ability to trap 
plasma, and this stuff is responsible for the observed 
edge-brightened features (mentioned in Section 1; Hard- 
castle et al. 2007). Relating to this, it is shown in Fig- 
ure 1 that the jet opening angle tends to decrease as 
away from the core; this trend continues up to knot G. 
This observational fact seems to suggest that such a mag- 
netic field certainly exerts the net collimation force for 
jet wide (Honda & Honda 2002), implying that the trans- 
verse correlation length of the magnetized filamentary 
turbulence reaches the radial size scale of the jet (beyond 
the aforementioned size scale of compact X-ray knots). 
Within the present framework, the various size scale of 
fragmented pieces inside the jet is considered to reflect 
the coherence length of the turbulence (as compatible 
with the argument given in Section 2.1), although the 
real situation might be more complicated (Blandford & 
Konigl 1979a,b; Hardcastle et al. 2003). 

Founded on this simple model, in the following Sec- 
tion 3 we attempt to expand rigorous arguments on the 
DSA of particles in the inner region (# pro j < 60"), in 
which the knots A and B have been, at least in part, 
associated with shock acceleration sites (Hardcastle et 
al. 2007). As the outer jet of # pro j > 60" concerns, the 
X-ray emission is largely diffuse, and the compact knots 
embedded in the diffuse emission arc much smaller than 
the radial sizes of the jet. Although there is no evidence 
for shocks (other than at knots) anywhere in the outer 
jet of # P roj > 60" at the moment, it seems that the ap- 
pearance of such inner structure as well as the measured 
X-ray spectrum can be explained by the present model 
including the inhomogeneity, according to the discussions 
given above. It is also noted that in the outer region, the 
presence of stochastic magnetic fields has been suggested 
(Mao & Wang 2007, within the framework of Fleishman 
(2006)). At this juncture, we provide, in the subsequent 
Section 4, a preliminary argument parallel to that ex- 
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pandcd in Section 3, considering an application of the 
model to the issue of particle acceleration in the outer 
region (60" < 6» proj < 200") including knots C, E, F, and 
G, where magnetic self-collimation of the jet seems to be 
viable (Figure 1). 

3. PARTICLE ACCELERATION IN THE INNER JET 
REGION 

3.1. Acceleration Timescale 

In the filamentary medium, the diffusion property of 
nuclei is quite different from that of electrons. For in- 
structive manner, let us begin with considering the cold 
limit, in which any species of charged particles are gyro- 
trapped in the (local) mean magnetic field permeating 
through a filament. For the DSA in this regime, the 
particles will suffer conventional gyro-resonant diffusion, 
going back and forth across the shock discontinuity to 
gain their energies. Then, the mean acceleration time 
is roughly proportional to <~ r g /c, where r g is the gyro- 
radius (smaller than the radial size of the filament) . In 
the ordinary case in which radiative mechanism sponta- 
neously impedes particle energization, the kinetic ener- 
gies disperse among the particle species, such that the 
ratio of the resulting gyro-radii among proton, arbitrary 
ion, and electron would indicate 



f 



A 



: r, 



9,e 



1: Vz(z 



(2) 



respectively, where Z and A are the charge and mass 
number of nuclei, and m e;P is the electron and proton 
mass, respectively. That is, we read off r 9iP ~ r Sj ,, and 
r g,p/ r g,e ~ O(10 6 ). The result indicates that the proton 
(ion) escape from magnetized filaments, which is asso- 
ciated with the bound-free transition (Honda & Honda 
2005), is easier than the electron escape. 

Unlike the major electrons quasi-secular ly bound to 
the filaments (Section 2.1; Honda & Honda 2007; Honda 
2008), energetic ions are likely to freely meander in the 
forest of filaments, undergoing the three-dimensional rms 
deflection. The latter is in the regime of the off-resonant 
scattering diffusion, which still makes possible to go back 
and forth across the shock, to additionally accelerate 
the transitionally injected ions to a higher energy range 
(Honda & Honda 2004, 2005). It is mentioned that this 
original idea could be reconciled with the recent simu- 
lation results for the similar plasma configuration (Sec- 
tion 2.1; Spitkovsky 2008). 

When considering the ideal Fermi type-I mechanism, 
the energy spectral index of accelerated ions ought to 
be identical to that of co-accelerated electrons (see, e.g., 
Schlickciser 2002, for details), and thus, denoted com- 
monly as p. The mean acceleration time for arbitrary 
nuclei is found to be given by 



\/Enai(0,p) c 



R 



E 



V 

ZeBU 



where a\ is the dimensionless factor 

/3(2p+l) 



ai(P,p) 



09-l)(p2-l)' 



(3) 



(4) 



which is of order unity, and e, E, R, and U are the 
elementary charge, particle energy, (transverse) correla- 



tion length of the filamentary turbulence, and flow speed 
upstream viewed in the shock rest frame, respectively 
(Honda & Honda 2005). Recalling the arguments given 
in Section 2.2, R is reasonably compared with the ra- 
dius of the jet, rather than that of compact knots. It is 
noted that Equation (3) is valid for the unbound parti- 
cles, which meander in the magnetized filaments, though 
still confined in the system with the size ~ R. Hereafter, 
we set p = 2.4, along with the arguments in Section 2.1. 
In this Section 3, we pay attention to the inner regions 
that wrap the well-resolved small X-ray knots AX1, AX2, 
AX3, AX4, BX2, and BX5 (Kraft et al. 2002; Kataoka 
& Stawarz 2005). For convenience, the parameter values 
(referred for the analysis in Section 3) are summarized 
in Table 1. Equation (3) is used to balance the various 
loss timescales evaluated below, to solve for E. 

3.2. Temporal Limits 
3.2.1. Shock- Accelerator Operation Time 

The dynamical timescales concerning the energy limit 
of accelerated particles include the operation time of 
shock accelerator and adiabatic expansion loss time. It is 
hard to precisely evaluate the former, because of the diffi- 
culty of the identification of shock structure and U . In an 
ideal case in which the shock is stationary as sustained by 
an exhaustless incoming flow, one can ignore this kind of 
restriction, so that the maximum energy analysis can be 
simplified (the related discussion is given later). Below, 
we consider, for heuristic, the generic situation in which 
the shock accelerator moves along the jet. For simplic- 
ity, it may be reasonable to compare U to the moving 
speed (with respect to the core). The shock accelerator 
operation time can then be estimated as t s h ~ L/U, and 
the adiabatic expansion loss time, i ac j ~ 3L/(2U r ) (e.g., 
Miicke et al. 2003), where L and U r are the length scale 
of the jet and the speed of radial expansion, respectively. 
For U > U r compatible with the self-collimating char- 
acteristic of jet flow (Section 2.2), we have the relation 
of t s h < tad, which suggests that i s h preferentially limits 
the particle acceleration. 

When taking into account of the relativistic effects, the 
scaling of t s h can be explicitly written as follows: 



t sh = 2.1 x 10 



10 



1 



L 0.5c 



r(l/)lkpc u 



(5) 



Here, T = [1 - (U / 'c) 2 ]~ 1/2 , and L is the deprojected 
distance of the concerned knotty regions from the core, 
that is, -L pro j/sin#, where L pro j and 9 are the projected 
distance and viewing angle, respectively (Section 2.1). 
According to the discussions given above, we set U — 
U app /[sin9 + (U app /c) cos 9], where c7 app stands for the 
apparent speed of motion. In Figure 2, we plot i s h as a 
function of 9, for the cases of J7 app = 0.1c (top) and 0.5c 
(bottom). It is noted that the enhancement of t s h for 
the marginally smaller 9 is due to the dominant increase 
of L against the increase of U and T in the right-hand 
side of Equation (5). For a plausible value of 6 = 15° 
(shaded lines; Hardcastle et al. 2003), one can find i s h ~ 
10 11 — 10 12 s. For example, for a possible parameter set of 
Uapp — 0.1c and 9 = 15°, we figure out £ s h = 3.5 x 10 11 s 
and 1.5 x 10 12 s in the regions with knots AX1 and BX5, 
respectively. 
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Viewing angle 0(deg) 



and U = U(U app ,0) (Section 3.2.1). 

Relating to this loss process, here it might be useful 
to remark the collision with photons, which also lim- 
its the particle energization. This actually takes place, 
as far as a finite it ra d level is sustained by the syn- 
chrotron radiation from co-accelerated electrons, galac- 
tic emissions, not to mention cosmic background lights. 
The timescale of the pr/ interaction can be estimated 
as t pl ~ (xf) ~ 1 (u m /u- ca _a)t syn , where x is a dimension- 
less factor almost independent of the source parameters 
(Biermann & Strittmatter 1987; Romero et al. 1996). 
By invoking \ ~ 200, thereby x^ ^ 0(1), we read off 
Vy ^ (u m /u I& d)t S yn, so that the inequality of u ra( j < u m , 
which is expected in the concerned jet-knot environments 
(Section 2.2), results in t pl > t syn . It thus appears that 
for the maximum energy analysis, it is sufficient to bal- 
ance f acc = i syn , as carried out above. 

In Figure 2 for U app = 0.1c (top) and 0.5c (bottom), 
we show the t syn levels for protons, i.e., Equation (6) 
for (A, Z) = (1,1), using the parameter values listed 
in Table 1. For smaller 9, t syn sharply increases, since 
6 increases (i.e., B decreases). As 9 increases, both 
6 and U decrease, so that their competition influences 
upon whether t syn decreases or not in the larger ^-region. 
For the concerned regions, the timescale ranges from 
t syn = 4.5 x 10 13 s (BX5) to 1.5 x 10 14 s (BX2) for 
9 = 15°. It is found that the t syn level is order of magni- 
tude higher than t sn over the whole #-range. As would be 
expected, the characteristic of the timescales well sepa- 
rated markedly simplifies the maximum energy analysis. 



Fig. 2. — Operation timescales of shock accelerator t s ^ for the 
apparent speed of Uapp = 0.1c (top) and 0.5c (bottom) vs. 9, 
in the regions around the six inner bright X-ray knots (included 
in the projection of pro j < 60"; labeled). The synchrotron 
cooling timescales t syn are also indicated (the explanatory notes 
correspond to those for t a h). 



3.2.2. Radiative Losses 

The radiative losses of the charged particles me- 
andering in the filamentary medium compete with 
the energization. We anticipate that the emis- 
sion is of diffuse (Fleishman 2006), and then, 
the cooling timescale that can be expressed as 
t syn = 36Tr 2 f(P)(A/Z) 4 [m p c 7 /(c 4 B 2 E)], where f(f3) = 

(6n)- 2 [f3(f3 + 2) 2 (/3 + 3)]/[2' 3 (/3 2 + 7/3 + 8)] (Honda & 
Honda 2007). For the maximum energy analysis, the 



3.2.3. Colhsional Loss 

The collision of the cosmic-ray (test) particles with tar- 
get particles also degrades the acceleration. In what 
follows, we clarify the condition for which the colh- 
sional loss can be neglected. The timescale of pp in- 
teraction may be simply given by t co \ = (a 



where a 



i'i> 



4 x 10" 



are the cross sec- 



tion and number density of target protons, respectively. 
In conjunction with the issue of neutrino production, 
the pp collision involves the possible reaction channel of 
p + p — ► 2c? + e + + v e + 7. Here we derive, for conve- 
nience, the critical density of target protons, denoted as 
n cr , above which t co \ < min(t s h, t syn ) ~ t sn sets in. 
In Figure 3, we plot n CI against 9. The n CT (9) depicted 



here is for the case of smaller U a 



0.1c, since the 



t syn = 8A x 10 



ai(P, P ) 1/3 
100 fiG 



value is lower (thereby, more crucial) than that for the 
t S yn gives the solution of E, larger U app case. For the inclination of 9 = 15°, the 

critical density ranges from 5.7 x 10 2 cm~ 3 (BX5) to 
n cr = 2.4 x 10 3 cm -3 (AX1). It turns out that, if n < 
10 2 cm -3 is satisfied, t co \ cannot be, in the whole #-range, 
the shortest timescale among the energy loss processes. 
For example, the expected value of the order of n 



temporal balance of t a 

and putting this into the i syn -expression denoted above 
is found to provide the following practical scaling 



m 



-1 2/3 



10- 



Z 2 / 3 \2Z J 



8/3 



B 



/100 pc 



1/3 



2/3 

-7- ) B. (6) 



0.5c 
U 



Note that setting to ((3,p) = (5/3,2.4), which reflects 
(a r ,a x ) = (0.7,1) (cf. Section 2.1), gives ai = 3.0 and 
f = 2.4 x 10~ 3 . In Equation (6), B is the quantity in the 
shock rest frame, and hence, for the moving shock case, 
set to B cq ^s = i5~ 5 / 7 (Kataoka & Stawarz 2005), where 
5 = r _1 [1 — (U/c) cos^p 1 (> 1) is the beaming factor, 



10 cm -3 (Meisenheimer et al. 2007) leads to t co \ > 
min(i s h,t S yn)' That is to say, one can almost ignore the 
collisional effects on the energy loss. It follows that the 
pp neutrino production channel is less outweigh than the 
P7 one that is considered in Section 3.5. 

Summarizing the discussions in Section 3.2, we can 
claim that the temporal balance equation can be reduced 
as i acc = min(i sh , i syn , t co \) — t sh , when considering the 
finite dynamical timescale of a moving shock. The solu- 
tion for E is denoted as E t , which indicates the maximum 
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Fig. 3. — Critical density n cr in each regions (labelled), above 
which the proton-proton collisional timescale t co i is shorter than 
*sh (f° r U app = 0.1c). For 8 = 15° (shaded line), it is evident that 
the collisional loss of an accelerated proton is less dominant, if 
n < n cr = 10 2 - 10 3 cm -3 . 



possible energy of an accelerated particle, determined by 
the temporal limits. This value should be compared with 
the energy restricted from spatial scale; the related issue 
is discussed below. 

3.3. Spatial Limit 

In the present context, when the mean free path of 
accelerated particles, transverse to the filaments (Aj_), 
reaches the order of the radial size of the jet (~ R), the 
particles would start escaping from the concerned sys- 
tem. This is just a subject to the spatial limit of particle 
acceleration. The balance equation X±(E) — R, con- 
trast to r g — R for the conventional gyro-resonant DSA, 
appropriately provides the solution for E (denoted as 
E s ), which indicates the maximum possible energy deter- 
mined by the spatial limit. Making use of the expression 
of Ax that can be derived from the quasi-linear kinetic 
equation (see Honda & Honda 2005, for the details), we 
obtain 



£ s =4.3x 10 18 a 2 (P,p)Z, 



D 



R 



100 100 pc 
where a 2 is the dimcnsionless factor 

-(/?-l)(p + 2)(p + 4)l 1/2 



eV, (7) 



a2(/3,p) 







(8) 



(5/3, 2.4), Equation (8) gives a 2 



Note that for (/3,p) 
3.4. 

If the outer ridge with the size exceeding R is per- 
vaded by a large-scale magnetic tail (Section 2.2), we 
could speculatively regard such a size as effective con- 
finement radius. When one allows for setting to , e.g., 
the value around ~ R ut(> R) listed in Table 1 (Kraft et 
al. 2002), the E s value (given in Equation (7)) is found 

to be raised by a factor of ~ ^jR out /R. However, this 
factor is of the order unity, and it can be confirmed that 
the effect does not severely affect the conclusions for the 



TABLE 2 

The Maximum Energies of Proton. Iron, and Neutrino 



Region 


£ p , m (EeV) 


£ 4>m (ZeV) 


£„, m (EeV) 


AX1 


35(130) 


0.92(3.4) 


1.8(6.5) 


AX2 


38 (130) 


0.99(3.3) 


1.9(6.4) 


AX3 


50 (140) 


1.3(3.8) 


2.5 (7.2) 


AX4 


43 (120) 


1.1(3.0) 


2.1(5.8) 


BX2 


64 (180) 


1.7(4.8) 


3.2(9.2) 


BX5 


120 (340) 


3.2(8.7) 


6.1 (17) 



Note. — Values in the round brackets indicate those for the 
stationary shock case. 



simple case excluding them. In this aspect, here we show 
the analysis, for which Equation (7) is referred in the 
rather conservative manner, to provide the theoretical 
upper limit of particle energy which is of the form that 
can be made a direct comparison with Et ■ 

3.4. Highest Energies of Proton and Iron 

For the concerned jet regions each, we evaluate the 
achievable maximum energy that is defined by E max = 
xjxm(E s , Et), where E 8 is given by Equation (7), and E t 
can be properly derived from i acc = < s h as argued in 
Section 3.2, to have the scaling of 



Et 



2.1 x 10 



L 



19 



1 



D 



ai(/3,p)V2' 



1 kpc 



1/2 



R 



100 pc 



100 fiG 

1/2 fjj^ 

0.5c 



1/2 



eV. (9) 



Hereafter, the solutions E ma , x for proton (A, Z) = (1, 1) 
and iron (56, 26) are, for convenience, simply denoted as 
E p<m and Ei, m , respectively. 

In Figure 4, for Z = 1 we plot E s (Equation (7)) and 
Et (Equation (9)), using the parameter values listed in 
Table 1. One can see that E PtTn = min(E s , E t ) = E t 
is ensured in the major range of 9 > 5°, while in 
the marginally smaller 9, the spatial limit regime with 
E p m — E a can appear. Particularly, for a standard 
6 = 15 c 



it is sure that E, 



E t is valid for both 
C^app = 0.1c and 0.5c in all concerned regions. This 
is one of the most noticeable results, which is different 
from the previous one derived from the analysis for a 
narrower jet with R < 100 pc (Vir A; Honda & Honda 
2004), where mm(E s , E t ) — E s appears in a possible U 
range. For practical reasons, the resulting values of E Pim 
for C/ app = 0.1c are listed in Table 2. It is noted that 
the Z7 app -dependence of E p ^ m is weak around 6 = 15°, as 
seen in Figure 4. We can conclude that the E p>m value 
more likely reaches as high as ~ 10 20 eV in the inner jet 
region (6 < 60") including knots BXn. Also, the values 
of E i m have been derived in parallel to E p m , and added 
in Table 2. It is found that the inner regions can be iron 
Zevatron; particularly, in the region with knot BX5, Ei >m 
reaches 3.2 x 10 21 eV. 

Here, we provide the discussion on the special case in 
which the shock is stationary (Hardcastle et al. 2003). 
The operation of DSA continues until the incoming flow 
into the shock breaks off; otherwise, the accelerator secu- 
larly operates. In the former case, it is hard to determine 
i s h, because of the complexity of the spatiotemporal dy- 
namical evolution of the jet. For the present maximum 
energy analysis, it might be adequate for considering the 
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Fig. 4. — Proton energies E a and Et determined by the spatial and temporal limits, respectively, in each regions (labeled). The axes are 
common in the panels. Note that Et(9) is the solution for min(i s i 1 , t syn , t co \) = t s h, given C/ a pp (labeled). The achievable maximum energy 
is given by E m ax = min(E s , Et); for both t/ a pp = 0.1c and 0.5c, E m ax = Et is satisfied at 9 = 15° (shaded lines), and for the stationary 
shock case, -Emax = E s . The -Emax values (denoted as E m ,p for proton) are summarized in Table 2. 



latter, and then, we read the temporal balance equation 
of ^acc = mm (iad> ^sym ^coi): which yields an enhanced 
value of E t , such that min(£' s , E t ) — E s more likely ap- 
pears. The -E s -levels evaluated from Equation (7), set as 
B = -E>oq,5=i, are indicated in Figure 4. For convenience, 
we list the resulting -E p m -values in Table 2 (in round 
brackets). In comparison with the results of t/ ap p = 0.1c, 
the values of E pnl increase, since the -E s -values are larger 
than E t (determined from £sh), as seen in Figure 4. 

The magnetic field strength involves the uncertainty, 
coming about through some assumptions in the esti- 
mation of the equipartition value (Section 2.2). Thus 
we also investigate the field intensity dependence of 
E p /i : m- In Figure 5, we show E PiTn against B eq< s=i for 
9 = 15°; the linear relations simply stem from E s / t oc 

B(= B cqi s=i5~ 5 / 7 ) in Equations (7) and (9). It is found 
that in the region with knot BX5, E Ptm (Ei )m ) > 10 20 eV 
can be achieved for the ranges of -B cq ,<5=i > 230 fiG(9 fiG) 
and > 85 /iG(3 ^G) for [/ app = 0.1c and stationary shock 
cases, respectively. 

3.5. High-energy Neutrino Emissions 

The p7 interaction discussed in Section 3.2.2 is revis- 
ited in light of the major neutrino production mecha- 
nism. The reaction sets in when the center of mass en- 
ergy of a j»7 interaction exceeds the threshold energy 
for the A-resonance; that is, E p , m is required for ex- 



ceeding the quantity of (m A — m p )c 4 / '4E y (= E Pyt h), 
where to a and E^ are the mass of the A-particle and tar- 
get photon energy, respectively (e.g., Halzen & Hooper 
2002). The threshold energy is found to scale as E Pyt h = 
1.6 x 10 17 r 2 (l cV/E-y) eV. Above E Ptth , charged' and 
neutral pions can be created by 737 — > A — > nn + and 
P7 — > A — > pir°, and the former triggers the decay chain 
of 7r + — > + — > e + + v e + Dy, + is to produce the 
very high energy neutrinos. 

Provided that the four final state leptons in the decay 
chain equally share the pion energy, the threshold neu- 
trino energy can be estimated as -E^.th = \ (xp^^) E P! th, 
where (xp^^) is the average fraction of energy trans- 
ferred from the initial proton to the produced pion 
(Halzen & Hooper 2002) . By virtue of the active galactic 
nucleus (AGN) unification scheme (cf. Section 5; Urry 
& Padovani 1995; Tsvetanov et al. 1998), it might be 
justified to adopt the estimation of (x p ^ w ) ~ 0.2, which 
is adequate for blazar jet environments. Then, for the 
target photon energy in the range of E 1 ~ 1 eV (for the 
CenA jet; e.g., Hardcastle et al. 2006), we work out at 
E^th — 8 x 10 15 eV. Similarly, one can estimate the 
achievable maximum neutrino energy (denoted as E Vt7n ) 
with reference to the E p _ m values obtained above. For 
convenience, the resultant E v _ m values are listed in Ta- 
ble 2. We find that E VyTa is likely to be of EeV range 
in the inner region; for J7 app = 0.1c and the stationary 
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Fig. 5. — Maximum proton energy E PtTn vs. the magnetic field 
strength compared to an cquipartition value B e „$- 1 (= B5 5 / 7 ), 
for U app = 0.1c and stationary shock cases. The explanatory 
notes for the former are the same as those given in Figures 2 and 3. 

case, the uppermost value is expected to reach as high as 
6 x 10 18 eV and 2 x 10 19 eV, respectively, in the region 
including knot BX5. 

Concerning the observability of the neutrino flux (fj) v ) 
by a /x-scintillation counter with the conversion proba- 
bility (P^p), the number of detected events could be 
expressed function of neutrino energy, such that 

-^events ~ 4>vPv^^ oc E v 1 (Halzen & Hooper 2002). 
Since we have little knowledge of the actual efficiency 
of energy conversion from accelerated protons to neutri- 
nos in FR-1 radio jets, the spectral normalization of <j) v 
involves the large uncertainty. Considering the Auger 
flux, however, Cuoco & Hannestad (2008) attempted to 
calculate the neutrino event rate in the detectors like 
IceCube, and found a rate of ~ 0.4 — 0.6 yr" 1 above a 
threshold of 10 14 eV. In any case, the estimated E vm 
provides the truncation energy of the neutrino spectrum, 
and in turn, could influence the neutrino flux estimation 
when the spectral index is given. The related issue is 
important for clarifying the share of cosmic-ray energy 
above the "knee" among celestial sources, but seems to 
be somewhat out of scope in this paper. 

4. ON THE PARTICLE ACCELERATION IN THE 
LARGE-SCALE JET 

4.1. Spatiotemporal Limits 

It may be worth applying the present DSA scenario to 
the unsolved issue of particle acceleration in the outer 
region more distant from the galactic core. Here, we 
preliminarily show the -E max distribution in the large- 
scale jet, simply providing B as a parameter (without 
reference to the values given in Table 1), and spatial 
profile of n. This might be more meaningful, if the length 



TABLE 3 

Summary of the Parameter and i? ma x Values 



Region 


L(kpc) a 


R(pc) b 




£ p , m (EeV) c 


£ Iim (ZeV) c 


A 


0.910 


352 


5/3 


6.5(25) 


0.17(0.66) 


B 


3.94 


819 


5/3 


21 (59) 


0.54(1.5) 


C 


4.93 


1310 


2 


33 (100) 


0.85 (2.7) 


E 


7.29 


1310 


2 


40 (100) 


1.0(2.7) 


F 


8.95 


2120 


2 


56 (170) 


1.5(4.4) 


G 


13.4 


2120 


2 


69 (170) 


1.8(4.4) 



a The deprojected length, assuming the viewing angle of 8 = 15° 
(Kraft et al. 2002; Hardcastle et al. 2003). 
b cf. Kraft et al. (2002). 

c Explanatory note for the round brackets is the same as that in 
Table 2. 



and radial scales of jet preferentially determine E maK , as 
we have seen this property in the inner region (Section 3). 
In the following, we investigate the condition, for which 
radiative and collisional losses arc negligible, to complete 
the discussions on the feasibility of UHECR production 
in the large-scale jet up to # pro j ~ 200" (£ pro j — 3.4 kpc, 
including the outer knots C, E, F, and G; e.g., Kraft et 
al. 2002; Hardcastle et al. 2007). 

For the analysis in Section 4, we adopt the parame- 
ter values listed in Table 3, with the use of the standard 
parameter set of l/ app = 0.1c and p = 2.4. We assume 
(3 = 5/3 in the inner region including knots A and B, 
and (3 = 2 in the outer region including knots C, E, F, 
and G, according to the arguments of the filamentary 
turbulent states (Section 2.1): in Equations (7) and (9), 
for instance, the latter case (/3,p) — (2,2.4), which re- 
flects (a r ,a K ) = (0.7,1.2), gives a\ — 2.4, a 2 = 3.8, and 
f = 4.3 x 10~ 3 . 

In Figure 6, we plot t sri and t syn , that is, Equa- 
tion (5) for given 9 as a parameter, and Equation (6) 
with (A, Z) = (1, 1) and (56,26), given -Beq,a=i as a pa- 
rameter (while 6 = 15° fixed). It is clearly seen that as 
L pro j increases, the value of i sn monotonically increases; 
especially for 9 = 15°, from t sh = 3.2 x 10 11 s (knot A) 
to 4.7 x 10 12 s (knot G). Note that the self-consistently 
determined t syn level for iron is comparable with that 
for proton, reflecting the physics that in the original 
timescale (being proportional to (A/Z^E -1 ), the en- 
hancement by the factor of (A/Z) 4 for iron is nearly 
canceled by the increase of E, which is caused by both 
the weaker synchrotron loss and stronger acceleration ef- 
ficiency (on account of the reduced i acc (oc 1/Z 2 )). Any- 
how, it appears that, as far as the possible condition of 
-Bcq,5=i < 500 /iG is retained, i s h is smaller than t syn in 
the entire region including knots A-G. 

As for the collisional loss, in parallel with the argument 
in Section 3.2.3, we provide the critical proton density 
n cr , above which t co \ < t s h sets in. For the timescale 
of nucleus-proton (Np) collision with t co \ ~ (a^pcn)^ 1 , 
we recall an empirical scaling of the cross section of 
OATp <~ 7rrQ^4 2 / 3 , where r w 1.4 x 10~ 13 cm. In Fig- 
ure 7, we explicitly depict n cr against £ pro j for pp and 
Np (say, A = 56) collisions, for the case of 9 = 15°. 
The level of n cr for Np is lower than that for pp, simply 
reflecting the thing that t co \ for Np is shorter than for 
pp (by the factor of <~ 0.04). The value of n cr decreases 
as iproj increases; e.g., for the pp case, it ranges from 
n cr ~ 3 x 10 3 cm" 3 (knot A) to 2 x 10 2 cm" 3 (knot G). 
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Fig. 6. — £ s h for given (labeled) and t Byn for given B ecl s=l 
(labeled) in each regions (labeled; up to the projected distance of 
~ 200" from the core). The open circles and diamonds indicate 
t S yn for proton (A, Z) = (1,1) and iron (56,26), respectively 
(providing 9 = 15°). Note that t sh < t Byn is satisfied for the con- 
ceivable ranges of 8 and B e q,<5=l throughout the concerned regions. 

Here, we attempt to make an ad hoc comparison with a 
density distribution deduced from the beta model that 
reasonably describes the galactic structure (Kraft et al. 
2002). By reasoning that the jet matter is concentrated 
on the opening angle of <f> (<~ 2R/L; i.e., the solid an- 
gle of 7r</> 2 /4), we invoke the crude density profile (along 
the jet) of n(Lp r oj) ~ (16/0 2 )n;,(r), where L WO} = rsmO 

and rib = no[l + ( r Ag) 2 ] _3( ^ 2 is the beta-model den- 
sity. In Figure 7, we exemplify the density profile of 
n '= I0 2 n b , reflecting 9 = 15° and (f> ~ 20° (Figure 1). 
Here, the parameter values of n = 4.0 x 10~ 2 cm~ 3 , 
tq = 0.5 kpc, and b = 0.4, have been adopted (Kraft et 
al. 2002). We then find that n takes the value in the range 
of <~ 10 _1 — 10° cm~ 3 in between the region of knot A to 
G (in accord with a value suggested by Meisenheimer et 
al. (2007)), and the level is order of magnitude lower than 
the n cr level for both the pp and Np collisions. Thus, it 
is claimed that the collisional effects on the energy limit 
are negligible. As a consequence, we can expect the sig- 
nificant reduction of min(i s h, t syn , t co \) <~ i s h, so that the 
temporal balance of t acc — t s h provides the correct solu- 
tion for E t . 

4.2. Highest Energies 

For each regions in the large-scale jet, we evaluate the 
achievable maximum energy of proton/iron, defined by 
E> p /i iT n = rnin(_E s , Et), where E s and E t are given by 
Equations (7) and (9), respectively. In Figure 8, we plot 
E s and E t (for Z = 1 and 26), providing 9 — 15° and 
B cq i5= i = 50 nG (Romero et al. 1996; Israel 1998; Hard- 




Distance from core L (arc min) 



Fig. 7. — Critical density n cr in each regions (labeled), above 
which t co i < i s h (for C/ a pp = 0.1c), for the proton— proton and 
iron— proton collision cases (labeled). For a comparison, the 
profile of the beta-model density nt, multiplied by factor 10 2 , is 
also depicted (short dashed) as a function of L pro j. For further 
explanation, see the text. 

castle et al. 2006; Meisenheimer et al. 2007) throughout 
the concerned regions. The assumption of the B e q,<5=i 
value, which corresponds to a spatially averaged one, is 
adequate for the present purpose, although the actual 
value might be larger, particularly, in the inner jet re- 
gion (Section 2.2). One can see that mm(E s , E t ) = E t is 
satisfied as a whole, for both proton and iron. The val- 
ues of E p n^ m are summarized in Table 3. It seems more 
likely that, in the outer jet region, E p>m can hardly reach 
~ 10 20 eV, even though setting i?-value to the uppermost 
range, whereas Ei^ m exceeds 10 20 eV, and reaches up to 
2 x 10 21 eV. When considering the stationary shock (in 
parallel to the arguments in Section 3), the E p /^ m value 
increases, and particularly, E p ^ m ~ 10 20 eV is achieved 
in the outer region. With regard to the neutrino pro- 
duction via the photopionization decay chains, we read 
off E v ^ m — j (xp^n) E PiTn > 10 18 eV in the outer region. 
Although the obtained results are preliminary, involv- 
ing parameter dependency, the expected overall features 
can be reconciled with the manifestation (given in Sec- 
tions 3.4 and 3.5) that the CenA jet is a candidate for 
the UHE proton accelerator that serves as EeV neutrino 
factory, and simultaneously, for the iron Zevatron. 

5. DISCUSSION AND CONCLUSIONS 

The validity of the filamentary model has been checked 
by reproducing the complicated blazar variabilities that 
involve the non- monotonous hysteresis patterns of broad- 
band electromagnetic spectrum in flare phases (Honda 
2008). By recalling the idea that regards FR-I radio 
sources as misaligned BLLac objects (Urry & Padovani 
1995; Tsvetanov et al. 1998), the application of the fila- 
mentary model to the CenA jet will be reasonably jus- 
tified. As for the particle acceleration in the filamen- 
tary medium, we point out that the present DSA mech- 
anism for nuclei, which owes to the (off-resonant) three- 
dimensional rms diffusion across the shock (Section 3.1), 
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Fig. 8. — Energies E s (open marks) and Et (filled marks) 
determined by the spatial and temporal limits, respectively, 
for proton and iron acceleration (labeled). Here, 6 = 15° and 
-B eqj( s=i = 50 fiG have been assumed throughout the concerned 
regions (labeled). The filled marks indicate the maximum energy 
obtained as -E m ax = mm(E s , Et) = Et for C/ a pp = 0.1c, where 
E s is indicated by the smaller marks; and the (larger) open 
ones indicate -Emax = E s for the stationary shock case. The 
Bmax values (denoted as E miP and £ m ,i for proton and iron, 
respectively) are summarized in Table 3. 

is simpler than the electron acceleration mechanism that 
is incorporated with the complexity of the energy hi- 
erarchy mediated by the transition injection (Honda & 
Honda 2007) . The acceleration of highest energy parti- 
cle in the CenA was first considered by Romero et al. 
(1996). In the early work, they dealt with the conven- 
tional gyro-resonant diffusion model even for proton, that 
is to say, implicitly supposed a large-scale ordered mag- 
netic field (with small- amplitude perturbations). They 
addressed that a proton could be accelerated to the en- 
ergy of 2.7 x 10 21 eV in situ, which is an order of magni- 
tude larger than the present E Pt , m values that are in the 
range of - 10 19 - 10 20 eV (Section 3.4; Table 2). Such 
a larger value was derived from simply equating a classi- 
cal acceleration timescale (e.g., Biermann & Strittmatter 
1987) with radiative loss timescales, but this fashion now 
appears to be somewhat optimistic (Section 3.2.2; Fig- 
ure 2). It can be claimed that taking into account of the 
shock operation time or particle escape is essential for 
the maximum energy analysis for the concerned source, 
to yield the reduced proton energy, which can hardly 
reach the 10 21 eV energy range. It is also mentioned 
that the situation in which the dynamical timescale lim- 
its the particle acceleration is analogous to the situation 
that typically appears in the supernova remnant shock 



acceleration, in which the adiabatic expansion of spher- 
ical shell likely becomes a major loss mechanism (e.g., 
Kobayakawa et al. 2002). 

In conclusion, we have evaluated the achievable maxi- 
mum energies of the proton and iron nucleus diffusively 
accelerated by the shock in the CenA jet including X-ray 
knots. In particular, we have taken into account of the 
more realistic DSA scenario that relies on the filamen- 
tary jet model responsible for the recent X-ray measure- 
ments, and elaborated the conceivable energy restriction 
stemming from spatiotemporal scales. The key finding 
is that, for the plausible ranges of the shock speed of 
~ 0.1c, viewing angle of > 10°, and magnetic intensity 
of < 500 fj,G, the uppermost particle energy tends to 
be inevitably limited by the shock accelerator operation 
timescale (< 10 13 s), rather than the radiative cooling 
losses (> 10 13 s). As a result, it has been demonstrated 
that there exists the acceleration region (off the galac- 
tic core), in which proton and iron can be energized to 
~ 10 20 eV and ~ 10 21 eV ranges, respectively. In partic- 
ular, for the inner region including the X-ray knot BX5, 
we work out at the maximum energies of 1 x 10 20 eV and 
3 x 10 21 eV for proton and iron, respectively; and esti- 
mate the corresponding cutoff energy of the neutrino pro- 
duced via the photopionization, as 6 x 10 18 eV. For the 
stationary shock case, we can expect the enhancement 
of the maximum energies, and read that the large-scale 
region up to the projection of about 200", as well, has 
a potential to energize proton and iron beyond 10 20 eV 
and 10 21 eV, respectively, to yield the pj neutrino with 
the energy exceeding 5 x 10 18 eV. We here manifest that 
the jet wide of the Cen A galaxy is a promising candidate 
for the UHE proton accelerator and Zevatron for high-Z 
nucleus, as desirable to account for the outcome from the 
Auger observatory (Abraham et al. 2007). 

The derived large values of particle energies are of those 
achievable at the acceleration sites. The observable en- 
ergies ought to be, of course, reduced, due to a major en- 
ergy loss mechanism including the photopion interaction 
with cosmic microwave background. However, the CenA 
is so nearby that the ballistic transport is anticipated to 
be not crucially degraded, on account of the weak de- 
cay property (such as — (l/E)(dE/dt) ~ 10~ 8 yr _1 for 
E > 10 20 eV; Romero et al. (1996)). We hope, in near 
future, the application of the present scenario to cos- 
mologically distant (super GZK) AGNs, to solidify the 
point source scenario for dominant UHECR production, 
which is responsible for the suppression of cosmic-ray 
flux above 4 x 10 19 eV that has recently been confirmed 
by the Auger's experiments (Abraham et al. 2008), and 
also, is of interest in conjunction with the detection of 
super GZK neutrinos (e.g., Ringwald 2006). 

I am grateful to Y. S. Honda for a useful discussion. 
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